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ABSTRACT 

We present simple models for disk evolution based on two different 
approaches: a forward approach based on predictions generic to hierarchical 
models for structure formation (e.g.. Mo, Mao, & White 1998) and a backwards 
approach based on detailed modeling of the Milky Way galaxy (e.g., Bouwens, 
Cayon, & Silk 1997). We normalize these models to local observations 
and predict high-redshift luminosities, sizes, circular velocities, and surface 
brightnesses. Both approaches yield somewhat similar predictions for size, 
surface brightness, and luminosity evolution though they clearly differ in the 
amount of number evolution. These predictions seem to be broadly consistent 
with the high-redshift observations of Simard et al. (1999), suggesting that 
the S-band surface brightness of disks has indeed evolved by ~ l.S'" from 
2; ~ to z ~ 1 similar to the models and is not an artifact of selection effects 
as previously claimed. We also find a lack of low surface brightness galaxies 
in several high redshift samples relative to model predictions based on local 
samples (de Jong & van der Kruit 1994; Mathewson, Ford, & Buchhorn 1992). 



Subject headings: galaxies: evolution 



1. Introduction 



Over the last several years, there has been a steady increase in the number and quality 
of observations available for disk galaxies from z = and z = 1. Schade et al. (1995,1996), 
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using early ground and space based images of galaxies from the Canada-France Redshift 
Survey (CFRS), found a net increase in the surface brightness of galaxies to 2; ~ 1. Along 
the same lines, Roche et al. (1998), compiling 347 galaxies from the Medium Deep Survey 
and other surveys, concluded that disk galaxies had undergone a net evolution in surface 
brightness and a net devolution in size. Lilly et al. (1998), using structural parameters 
extracted from HST images of the combined CFRS and LDSS2 sample, concluded that 
there has been essentially no evolution in large disks out to z ~ 1. As a preliminary effort as 
part of the DEEP survey, Vogt et al. (1996,1997) found little evolution in the TuUy-Fisher 
relationship (< 0.3™) out to 2; ~ 1. More recently, these observations have been augmented 
by the DEEP sample with 197 galaxies from the Groth strip to / < 23.5, 1.5 magnitudes 
deeper than the LDSS2-CFRS sample. In a first paper, Simard et al. (1998) concluded that 
there had been little evolution in the disk surface brightness distribution to 2 ~ 1 contrary 
to previous claims. 

A number of different approaches have been proposed for making specific predictions 
about disk evolution. Mo, Mao, & White (1998a) showed how the standard paradigm 
for hierarchical growth of structure combined with simple assumptions about angular 
momentum conservation led to simple scaling relationships for the change in disk properties 
as a function of redshift. Other authors (Ferrini et al. 1994; Prantzos & Aubert 1995; 
Prantzos & Silk 1998; Boissier & Prantzos 1999; Chiappini, Matteucci & Gratton 1997), 
taking more of a backwards approach to the problem, used detailed studies of the profiles 
of the Milky Way and other nearby galaxies to propose radially dependent models of star 
formation in disk galaxies, models which could be used to make detailed predictions about 
high-redshift disk evolution. 

Already there have been a number of elegant studies in which both the backwards 
approach (Cayon, Silk, & Chariot 1996; Bouwens, Cayon, & Silk 1997; Roche et al. 1998) 
and the forwards approach (Mao, Mo, & White 1998; Steinmetz & Navarro 1999; Contardo, 
Steinmetz, & Fritze-von Alvensleben 1998; van den Bosch 1998; Mo, Mao, & White 1998b) 
have been used to interpret the observations available for disk galaxies, mostly to 2; ~ 1. 

Unfortunately, none of these studies considered the important effect that a large spread 
in surface brightness could have on the interpretation of these observations, particularly 
the potentially large fraction of low surface brightness galaxies. In some studies, the 
surface brightness selection effects at low and high redshift were simply ignored, and 
in others, e.g., Roche et al. (1998), the spread was limited to 0.3mag/arcsec^ about 
Freeman's law (Freeman 1970). Clearly, given the apparent large numbers of low surface 
brightness galaxies seen locally, it is quite logical to wonder if these galaxies are detectable 
in current high redshift surveys. Indeed, one might wonder whether these galaxies or the 
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observed correlation between luminosity and surface brightness may have already affected 
the interpretation of high redshift observations. In light of the recent claim by Simard 
et al. (1999) that the apparent surface brightness evolution thus far inferred to 2 ~ 1 
is completely due to surface brightness selection effects, such a study would seem to be 
especially timely. Secondly, none of these studies directly compared the predictions of 
the forward and backward approaches using the same observations. Simple comparisons 
of the scaling expected in surface brightness, size, luminosity, and number are useful for 
interpreting the high redshift observations. 

To address these shortcomings, we shall therefore consider implementations of both 
approaches, normalize them to the observed 2 ~ size-luminosity relationship, compare 
their predictions, and consider how each of them fares at explaining the observed disk 
evolution out to z ~ 1 incorporating all the selection effects as they are best understood. 
We commence by presenting our models (§2) and the observational samples with which 
we compare (§3). We present the results (§4), discuss them (§5), and then summarize our 
conclusions (§6). Throughout this study, we use Hq = 50km/s/Mpc unless otherwise noted. 



2. Models 

We begin by sketching the base z = model to be used for both the models which 
follow. We use a set of gaussian LFs based on those presented in Binggeli, Sandage, and 
Tammann (1988): 

0(M)rfM = expi-i ^ ~^* f)dM (1) 

2TTaM o-M 

We adjusted the bulge-to-total (B/T) distributions of these galaxy types to obtain fair 
agreement with the de Jong & van der Kruit (1994) sample. Finally, we adjusted the 
luminosity function so there was rough agreement with the combined Sabc and Sdm 
luminosity functions presented in Pozzetti, Bruzual, & Zamorani (1996). We present our 
parameterized populations in Table 1. 

For the above luminosity functions, we convert z = i?-band luminosity to mass using 
a constant mass-to-light ratio, where the mass of a Mf^j = —21.1 galaxy is 1.1 ■ IO^^Mq. 
We assume a log-normal scatter of 0.3 dex to reproduce the observed Tully-Fisher scatter 
though variation in the formation times (van den Bosch 1998) and concentration indexes 
(Avila-Reese, Firmani, & Hernandez 1998) certainly play a role. 

To translate this mass into a circular velocity and size, we calculate the time at which 
the ambient halo formed. Since halos are always accreting more mass and merging with 
larger halos, there is some ambiguity in defining this, so for simplicity we take it to equal 
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the redshift at which half of the mass in a halo has been assembled. We determine the 
distribution of formation times using the procedure outlined in Section 2.5.2 of Lacey & 
Cole (1993). We take the circular velocity of the halo to be that corresponding to the halo 
at its formation time using Eq. (14). 

Then, given the mass, circular velocity, and luminosity of the z = disk, we randomly 
draw the sizes rg from the following distribution: 

(JaVStt 2 crA/ln(10) 

where r* = 6.9kpc, ax = 0.37, and M^^g = —21.1. We provide a basic observational and 
theoretical motivation for this scaling in §2.1. Note that here the surface brightness is 
proportional to L^^'^, that the spread in the size distribution is proportional to the spread 
in the distribution of A, and that the scale length (surface brightness) of the average L^, 
galaxy is exactly equal to that predicted by Freeman's law. 

We assume that the SED of disks and bulges is identical to that of a 10 Gyr-old 
stellar population with an e-folding time of 4.5 Gyr and a. tb = 0.3 foreground dust screen, 
the extinction curve being that of Calzetti (1997). For simplicity, we assume this SED is 
constant independent of time. We use the Bruzual & Chariot tables from the Leitherer 
et al. (1996) compilation for this calculation. In the following models, we evolve the size, 
number, and luminosity of all galaxy types using simple single- valued functions of redshift: 

R{z) = R{0)Er{z) (3) 

N{z) = N{0)Er,{z) (4) 
L{z) = L{0)El{z) (5) 



Table 1. Model parameters for disk population. 



(po (^M B/T aB/T 

2.0 1.1 -19.4 0.25 1.8 

2.5 1.1 -19.6 0.08 0.5 

8.0 1.3 -18.4 0.04 0.25 

24.0 1.3 -16.0 0.01 0.25 
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For simplicity, we assume similar scalings in the properties of bulges as a function of time. 

Since the color and luminosity of disks has been shown to be correlated with inclination, 
it is reasonable to suppose that disks are not transparent. Unfortunately, there is much 
controversy concerning the degree to which disks are or are not transparent. For better 
or worse, we will side-step this controversy and simply adopt the Tully & Fouque (1985) 
prescription for extinction in the B band: 

An = -2.51og(/(l - exp(-rsecz)) + (1 - 2/)(l^^^^PizI^)) (g) 

r sec I 

where r = 0.55, / = 0.25, and i is the inclination of the disk, corresponding to a face-on 
disk. We shall assume our model galaxies are always observed at an inclination of 70 deg 
and therefore always correct the observations to this inclination for comparison with the 
models. In the B band, this corresponds to an extinction correction of 0.67"*. 



2.1. Hierarchical Model (Forwards Approach) 

The use of simple scaling relationships between the properties of disks and the halos 
in which they live has provided a relatively successful way of explaining both the internal 
correlations between disk properties and their evolution to high redshift. In this picture 
developed by Fall & Estathiou (1980) and revived more recently by Dalcanton et al. (1997) 
and Mo et al. (1998) among others, the bivariate mass and angular momentum distribution 
nicely translates into a luminosity and surface brightness relationship for disk galaxies, mass 
translating directly into luminosity and the dimensionless angular momentum translating 
directly into surface brightness. 

This picture has had much success in explaining the internal correlations between size, 
circular velocity, and mass. For example, De Jong & Lacey (1999) recently showed that 
the observed local bivariate luminosity-size distribution is nicely fit by this picture, albeit 
with a slightly smaller scatter in surface brightness than might otherwise be expected. 
For constant mass-to-light ratios, the rough M relationship for halos provides a 

relatively natural explanation for the luminosity-circular velocity (TuUy-Fisher) relationship 
(Dalcanton, Spergel, & Summers 1997a; Mo et al. 1998a; Steinmetz & Navarro 1999; 
Contardo et al. 1998). Finally, the R oc Vc relationship found in galaxy halos is similarly 
observed in the disk population (Courteau 1997). 

On the other hand, this picture says little, if anything, about how the gas disk 
evolved over time and therefore what its local properties (i.e., spatial variations in the 
metallicity, color, stellar ages, gas density, etc.) or global properties (total gas mass) are. 
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and therefore comparisons of this sort will depend upon the model adopted, whether it 
be one of the popular semi-analytic approaches (Cole et al. 1999, Somerville & Primack 
1998) or a full N-body hydrodynamical simulation (e.g., Contardo et al. 1998). Moreover, 
it is now apparent that the simple scaling model is fundamentally flawed with regard to 
the implementation of galaxy formation theory as revealed by high resolution numerical 
simulations (Moore et al. 1999; Navarro & Steinmetz 1999; Steinmetz & Navarro 1999). 
Since it however is the only detailed model available, it is imperative to fully explore 
comparisons with data, properly incorporating observational selection effects, in order to 
establish the correct basis for ultimately refining the model. 

For a detailed discussion of this picture, i.e., the idea that simple scalings in the 
properties of halos lead to simple scalings in the properties of disks, the reader is referred 
to Mo et al. (1998a) and later papers by the same authors. For the sake of clarity, we shall 
review some of this material. In the standard spherical collapse model for an Einstein-de 
Sitter universe, the density of the collapsed halo is 187r^ ^178 times the critical density 
of the universe at collapse time (see also Gunn & Gott 1972; Bertschinger 1985; Cole & 
Lacey 1996), but depends on the density of universe through the parameter x = 1 — Q{z). 
Expressing the result in terms of the super critical density parameter 

T-^ = ^cPc (7) 

3 " ' vir 

Bryan & Norman (1998) found that for Q + Q\ = 1, 

Ac ^ ISn^ + 82x - 39x2 ^^g^, 

for fl^ = 0, 

Ac ^ I8n^ + 60a; - 32x^ (9) 

where x = Q{z) — 1. 

Using the virial theorem, it is possible to write equations to relate the mass, radius, 
and circular velocity of each halo. As in Somerville & Primack (1998), it can be shown that 

Kl =™- ^H{z)\l^ (10) 

where r^i^ is the halo size, Vmr is the circular velocity of the halo at r^j^, G is Newton's 
constant, and 



H{Zf) = //oV^A.O + (1 - f^O - f^A,0)(l + Zff + + Zff. (11) 

Using the fact that M = AcPcl^rr^j^ = AcH{z)^^rf,^^, we can rewrite this as 

Kl = l{A,-n^)H{zrrl, (12) 
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or 

'^"'^ (13) 

/l(A,-fiA)i7(2) 

Similarly, we can now rewrite the halo mass as 

1/2. r ■ V^- 

G GH{z)^\{A,-n^) 

Assuming the matter which settles in the disk to be some fraction of the halo mass and 
the angular momentum of this settling matter to be some fraction jd of the halo's angular 
momentum, a straightforward derivation (e.g.. Mo et al. 1998a) allows one to obtain 

M, = ""f^"- (15) 

for the mass of the disk and ^ 

Rd = —f^ ( — ) Ar„ir (16) 

for the radius of the disk. The dimensionless angular momentum parameter A is defined as 

A = J\E\^''^G~^M-^''^ (17) 

where J is angular momentum, M is the mass, and E is the total energy of the bound 
system. 

There are three relatively simple reasons to go beyond this simple approach. First, the 
adiabatic contraction of the halo due to dissipation of baryons towards the halo center will 
modify the halo profile. Second, numerical simulations show that model halos actually have 
a Navarro, Frenk, & White (1997) profile for the the dark halo rather than the isothermal 
profile used above. Finally, in order to make comparisons back to the observations, it is 
important to consider the observationally-measured rotational velocities of the disk rather 
than the rotational velocities of the halo proper. Mo et al. (1998a) have found approximate 
fitting formulas for the consequent corrections made to the disk radius Rd and the circular 
velocity at 3/?^: 

i?, = i= (^) \r,„f~'/'fn (18) 

K(3i?,) = V,„fv (19) 
where approximate fitting functions for /c, /r, and /y are given by 

^'4-^ (215)°' P«) 
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/ ^ \ -0.06+2.71md+0.0047/A 

— j (l-3md + 5.2m^)(l- 0.019c + 0.00025c2 + 0.52/c) (21) 

. . -2.67^,-o.oo38/A+o.2A ^ 2^ 1 + 0.057c - 0.00034c2 - 1.54/c , , 

/v^ 77T 1 + 4.35m, -3.76m^ --— 22 

\0-l/ [-c/(l + c) +ln(l + c)]'/^ 

where c is the standard halo concentration parameter for the Navarro et al. (1997) profile. 
For simplicity, we use = 0.05 and c = 10 to convert Kir to Vc in order to compare with 
the observations. We do not use Eqs. (15)-(16) for these comparisons. Note that larger 
values of render disks unstable at relatively faint surface brightnesses and thus have 
difficulty accounting for Freeman Law-type surface brightnesses (Freeman 1970). 

On the basis of these simple halo scaling relations, the size and mass of disks at any 



redshift simply scales as l/iJ(z/)y |Ac(z/) — A(zj), being the redshift at which these 
high-redshift disks formed. Consequently, the size and luminosity scale as 

r(z) = — (23) 

L{z) = ^Az)M{z) = ^iz) , ^^'^ fv (24) 

^ ^ 'GH{zf{z))^\{AAzf{z))-n.{zs{z)) 

Using Eq. (3) and (4), we now have our functions Er{z) and El{z): 

^ , , ^(^/(0))v/|(Ac(^/(0))-fiA(./(0))) 

Er{z) = ' , = (25) 

H{zf{z))^l{AAzj{z))-Qj,{zj{z))) 

7(0)i7(^K0))v/|(Ac(^/(0))-fiA(^/(0))) 
= ,^ = (26) 

^{z)H{zf{z))^'^{AAzf{z)) - n^izfiz))) 

where 7(2;) is the mass-to-light ratio at redshift z. Note that this is quite different from 
scaling these disks simply in terms of the redshift at which these disks were observed, 
particularly in the case of low Q where little evolution in the size or baryonic mass of the 
disk population is expected. 

We assume that the z = luminosity function scales in number as a function of z in an 
analogous way to how the lO^^M© halos scale in number. Using the Press-Schechter (Press 
& Schechter 1974) mass function 

iV(M, z)dM = -x^^.^J^ exp (- ^ , \ ^^dM (27) 
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we see that the halo number density scales as 

noc[D{z)]-^ (28) 

since the exponential factor remains approximately unity for the IO^^Mq mass scale. In 
terms of the formalism of Eqs. (3-5), 

E.(.) = ^ (29) 

Here, D{z), the growth factor, was computed using the formula tabulated in Carroll, Press, 
& Turner (1992). 

We now provide a theoretical and observational justification for our size-luminosity 
distribution. Theoretically, in the Fall & Estathiou (1980) picture, the spread in surface 
brightnesses derives from the spread in dimensionless angular momenta for halos. An 
approximate parameterization of the dimensionless angular momentum distribution is 



ln(A/A)- 



For A = 0.05 and ax = 0.5, the above expression closely approximates the distribution 
obtained from N-body simulations (Warren et al. 1992; Cole & Lacey 1996; Catelan & 
Theuns 1996) and analytical treatments (Steinmetz & Bartelmann 1995). The above spread 
in dimensionless angular momentum directly translates into the following distribution of 
sizes: 

In (r/rg 



(j)(r)dr = —j= exp 



2al 



dlogr (31) 



For disks with a constant mass-to-light ratio, it follows from Eqs. (13-16) that 

"Td oc Tyir oc VmrH{zf)~^ oc M^^ H{zf)~^^^ oc L]/^ H {z f)'"^^^ . Ignoring the dependence of 

H{zf)~'^^^ on the luminosity, it follows that the surface brightness {Ld/r'^) scales as L]/^- 

In fact, de Jong & Lacey (1998) found that the Mathewson, Ford, & Buchhorn (1992) 
data set gave a good fit to the following bivariate size-luminosity distribution with similar 
properties to those predicted above: 

^/ M *o , l.logre/r! -0.4(M- M,)(2//3- l),2^ 
^r,,M)dlogrJM = L= exp(--[^^^^ ) '-f) 

where $o = 0.0033Mpc-^ a = -1.04, = 3, M, = -22.8, r* = 7.9kpc, and ax = 0.37 
(converting their sizes and luminosities from ho = 0.65 to the Hq = 0.50 used here). Implicit 
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in the above bivariate distribution is a distribution in sizes analogous to Eq. (31), a 
Schechter distribution in luminosity, and a SB oc L^^'^ correlation between luminosity and 
surface brightness. Similar scalings are apparent in the McGaugh & de Blok (1997) sample. 

Now let us compare a typical galaxy in this model with the observations. Using 
our stated assumption, a galaxy has a mass of 1.1 • lO^^M©. A typical formation time 
occurs aX z = 0.3. Using Eqs. (13-14), the circular velocity and size of the halo is 132 km/s 
and 272 kpc (compared to the 140 km/s and 241 kpc predicted assuming a constant 200pc 
for the collapse density as in Mo et al. 1998). Then, using Eq. (18), the size of the disk is 
~ 6.0 kpc. 

This is smaller than the empirical findings of de Jong & Lacey (1998) (7.9 kpc), our 
own comparisons to local observations (§4.4) (6.9 kpc), and Freeman's Law, which gives 6.9 
kpc. Supposing this to be due to a slight cut-off at low values of the dimensionless angular 
momentum due to disk instabilities (Efstathiou, Lake & Negroponte 1982; Dalcanton et al. 
1997; Mo et al. 1998a; van den Bosch 1998), we scale up the size of a typical L^, disk galaxy 
to 6.9 kpc and reduce the spread in dimensionless angular momenta to ax = 0.37, as found 
by de Jong & Lacey (1999) in the analysis of the Mathewson et al. (1992) sample. 

Throughout our analysis, we shall take the Q = 0.3, Qa = 0.7 model as our preferred 
fiducial hierarchical model because of its better correspondence with the evolution in the 
number of small disks observed up to 2; ~ 1 (Mao et al. 1998). We evolve the mass-to-light 
ratio 7(2;) as (1 + z)~°'^ to reproduce the observed evolution in the TuUy-Fisher relationship 
(see §4.2). 

2.2. Infall Model (Backwards Approach) 

Instead of trying to determine how the global structural properties of disks evolve 
based on the corresponding properties of their ambient halos, it is also possible to 
examine a number of local disk galaxies in great detail and to use detailed models of their 
observed properties (gas profiles, stellar profiles, metallicity profiles, current SFR profiles, 
age-metallicity relationships) to determine how galaxies might have evolved to high redshift. 

There is no consideration of how individual halos might evolve backwards in time in 
these models, both for simplicity and because of large uncertainties in the local distribution 
of dark matter. Consequently, while for the forwards approach, the entire evolution of disk 
properties derives from an evolution of the halo properties, the infall models considered 
here completely ignore these effects. Conversely, while the forwards approach presented 
here ignores issues related to the manner in which halo gas is converted into stars, for the 
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infall model, such issues are important. 

Naturally, given that one always adopts the observed local universe in this approach 
as known, this approach does not explain in and of itself why the local disk population is 
as it is. Indeed, it cannot since there is no link to the initial conditions. In this view, for 
the hierarchical approach, we adopted the local universe we did because it was a natural 
prediction of the model, and for the infall approach, we adopted it because it agrees with 
the observations. 

We examine such a model for the evolution of local galaxies based upon the Prantzos 
& Aubert (1995) model for the star formation rates, metallicites, stars, and gas content 
for the Milky Way disk. We previously presented this infall model elsewhere (Bouwens 
et al. 1997; Cayon et al. 1996), and we shall revisit it here. This model ignores radial 
inflows for simplicity and takes the star formation rate to be proportional to both the gas 
surface density (S^) and the reciprocal of the radius r, which for a fiat rotation curve is 
proportional to the epicyclic frequency: 

dt - r,{r) ^^'^ 

where Tg{r) = [0.3 (r/r©)"^ Gyr]~^. Physically, such a star formation rate results if the star 
formation rate is proportional to the rate at which molecular clouds collide (Wang & Silk 
1994) or the periodic compression rate (Wyse & Silk 1989). 

For simplicity, the accretion time scale Tff was taken to be independent of radius since 
a variation in this time scale is not strongly constrained by the observations (Prantzos & 
Aubert 1995). The spread in Hubble types was then naturally taken to arise from a spread 
in this time scale (Cayon et al. 1996). The equation for the evolution of the gas density is 
then 



d^g{r) _ ^^(r, T) + S,(r, T) e"*/^// S 



(34) 



dt \-e~^hf Tff Tg 

where T is the time from the formation of the disk to the present. Integrating these 
equations yields the result 



g 



Given the fact that this model derived from only one galaxy, it is difficult to know how 
to extend its evolutionary predictions to galaxies with different luminosities and surface 
brightnesses. One possible means of extending this model to galaxies beyond the Milky Way 
involves simply scaling the star formation rates by the differential rotation rate. The change 
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in the star formation rate would then be proportional to Vd R or the typical dynamical 
time for the disk. Since R is roughly proportional to V^, there would be no large change 
in the time scales as a function of disk mass or rotational velocity. Accordingly, Bouwens 
et al. (1997) simply elected to scale everything in size to reproduce all luminosities while 
conserving surface brightness, which is precisely what we have done here. 

To determine the scaling relations for galaxies using the infall approach, we performed 
the calculation for each galaxy on a series of 30 different rings varying logarithmically in 
size, where the smallest is a circle with radius 0.2 kpc and the largest is a ring of radius 
60 kpc with width 12 kpc as done in Bouwens et al. (1997) and Roche et al. (1998). We 
calculate the evolution in rest-frame B band magnitudes by evolving each ring separately 
to keep track of its gas mass, stellar composition, and metallicity, and we output its colors 
using the Bruzual & Chariot instantaneous-burst metallicity-dependent spectral synthesis 
tables compiled in Leitherer et al. (1996). 

In the rest-frame B band, we found the following scaling relationships: 



Clearly, without number evolution, we take En{z) = 1. 

For the sake of clarity, we note that the present model differs from the one presented in 
Bouwens et al. (1997) in terms of both the luminosity functions used and the bulge-to-total 
distribution assumed. Furthermore, in the Bouwens et al. (1997) study, the preferred values 
of the age T and gas-infall time scale r// used in the Prantzos & Aubert (1995) study 
were scaled to reproduce the number counts. No such scaling was attempted in the present 
model and we simply use the same Tff for all disk types. 

This model is similar in spirit to the size-luminosity evolution model presented by 
Roche et al. (1998) based on the infall models of Chiappini et al. (1997), which models the 
infall, star formation, and chemical evolution of both the thin and thick disk components. 
For the purposes of illustration, we shall compare the Chiappini et al. (1997) infall model to 
the one just described in §4.1, after which we will restrict our consideration to the redshift 
scalings given by the Prantzos & Aubert prescription. In this model, the star formation 
time scale is equal to 



The star formation commences at t/orm = (16Gyr — 0.35r). Note that the time scale for 
star formation here depends on the radius to the first power as in our model, the preferred 
Prantzos & Aubert (1995) model, and the recent work by Boissier & Prantzos (1999). 



El{z) = 10 



■0.4(-0.6z) 



(36) 



Er{z) = 1 - 0.27z 



(37) 




(38) 
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3. Observations 

3.1. Low-Redshift Samples 

We make comparisons against local samples with information on the size, luminosity, 
and circular velocity of local galaxies. Though in principle we could have just used the 
Courteau (1997) sample, we follow Mao et al. (1998) in using a compilation of three 
different samples for the comparisons which follow to examine the three two-dimensional 
relationships. 

de Jong & van der Kruit Sample: 

The de Jong & van der Kruit (1994) sample provides a nice sample for examining the 
local size/magnitude relationship. It is selected from the Uppsala Catalogue of Galaxies 
(Nilson 1973, hereinafter UGC), over only ~ 12.5% of the sky and uses only relatively 
face-on {h > 0.625) galaxies (37.5% of all orientations). Following de Jong (1996), we also 
take it to be diameter-limited in R to galaxies larger than 2' at 24.7 i?-band mag/arcsec^. 
Whereas de Jong & van der Kruit (1994) in their treatment of their sample assume 
transparent disks, we correct observed magnitudes to an inclination of 70 deg using the 
TuUy & Fouque (1985) inclination corrections. 

Courteau (1997) Sample: 

We use the Courteau (1997) sample to calibrate the local z = Vc — size relationship. 
The Courteau (1997) sample contains 304 Sb-Sc galaxies from the UGC with Zwicky 
magnitudes tub < 14.5, i?-band angular diameters larger than 1', and i?-band major axis 
< 4'. We take their Vopt = V"c(3.2i?d) as the circular velocity and the 25 r mag/arcsec^ 
isophote as the radius. 

Pierce & Tully (1988) Sample: 

We use the Pierce & Tully (1988) sample to calibrate the local z = Q — luminosity 
relationship. The Pierce & Tully sample was taken from galaxies in the area of the Ursa 
Major cluster and is complete up to Bt < 13.3. It includes all galaxies which are not 
elliptical or SO, not more face on than 30 deg, and not possessing confused H I profiles. Note 
that in this study and in the Vogt et al. (1996,1997) studies to be discussed, the observed 
absolute -B-band magnitudes were corrected to intrinsic (unextincted) values using the 
Tully & Fouque (1985) inclination corrections. 
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3.2. High-Redshift Samples 



Simard et al. (1999) Sample: 

For the magnitude-radius relationship, we use the data presented by Simard et al. 
(1999). This data set contains structural information for ~ 200 galaxies to / < 23.5 from 
6 different WFPC2 pointings in the Groth strip (~ SOarcmin^). Spectra were obtained for 
only a fraction of the faint galaxies, but for galaxies with spectra, there was nearly 100% 
redshift identification. Following Simard et al. (1999), we can quantify the selection effects 
of this sample. The probability that a galaxy with apparent magnitude /814 and radius 
would fall in the photometric sample is Sup{hi'i) fd)- From Figure 4 of Simard et al. (1999), 
we have approximated this as 



(Actually, this provides a steeper surface brightness cut-off than the selection function 
given in Simard et al. 1999.) The probability that a galaxy with apparent magnitude /814 
and radius would be selected from the photometric sample for spectroscopic follow-up is 
given by >S'p5'(/8i4, r^), which we have approximated as 



(again by eyeballing Figure 4 of Simard et al. 1999). Putting these two selection effects 
together, the probability of selecting a galaxy with apparent magnitude J814 and radius 
Td is simply the product of these quantities, namely, Sup{hi4:^fd)Sps{hi4:^fd)- Given our 
ignorance about the inclinations used in the Simard et al. (1999) study, we assume an 
average inclination of 60 deg in transforming the absolute magnitudes to an inclination of 
70deg using the Tully & Fouque (1985) law, which yields a correction of 0.27™ here. 

Lilly et al. (1998) Sample: 

We also use the data from the LDSS2-CFRS sample with HST WFPC2 follow-up 
to look at the magnitude-radius relationship. From Table 3 of Brinchmann et al. (1998), 
the effective of the CFRS portion of this sample is 0.01377 deg^ (49 arcmin ). The 
survey is magnitude limited to 17.5 < I < 22.5, where the magnitudes are isophotal to 
28.0 Iab mag/arcsec^. Though the surface brightness limit is quoted as 24.5 Iab mag/arcsec^, 
we have used the more conservative surface brightness limit 23.5 /^b mag/arcsec^. In 



' 1, /8i4 + 51ogrd(//) < 21, 

< 1- f(/8i4 + 5 log rd (//)), 21< J8i4 + 51ogrd(//) < 22.5, 
^ 0, 22.5 < /8i4 + 51ogrd(//). 




814 ~ 



hi4 < 19.3, 
19.3), 19.3 < J814 < 23.5 
23.5 < ^814 
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principle, then, for any galaxy detected the isophotal magnitude should be approximately 
equal to the total magnitude. Lilly et al. (1998) chose to examine that subset of galaxies 
from this data-set which were disk-dominated and had disk scale lengths > 4 kpc, a sample 
we shall henceforth refer to as the Lilly et al. (1998) large disk sample. Given that the 
central surface brightness is not strongly correlated with inclination angle, Lilly et al. (1998) 
concludes that disks are consistent with being opaque. We use the TuUy & Fouque (1985) 
to transform the listed absolute magnitudes to an inclination of 70 deg. 

Vogt et al. (1997) Sample: 

In contrast to high-redshift samples with both magnitude and radial information, 
high-redshift samples with circular velocity measurements are considerably smaller and 
possess less well-defined selection criteria. In fact, the Vogt et al. (1996,1997) sample with 
16 galaxies is the largest such published sample. The selection criteria for this sample is 
still somewhat qualitative and patchy in nature. It considers galaxies with an inclination 
greater than 30 deg, detectable line emission, undistorted disk morphology, and an extended 
profile. We assume an /814 < 22.5 magnitude limit as used in the Vogt et al. (1997) sample. 

4. Results 

4.1. Basic Scalings 

Before getting into a detailed comparison of the models with the observations, we 
begin by illustrating the manner in which the sizes of (M ~ 1.2 • lO^^Mo) galaxies 
typically evolve as a function of redshift for the different models in Figure 1. We present 
this evolution in terms of the rest-frame B half-light radius, the specification of a band and 
a measure being necessary to the intrinsic band and profile dependence of evolution in the 
infall model. 

The hierarchical models predict more size evolution than expected from both infall 
models considered here. Of course, the f2 = 1 model possesses more size evolution than 
the VL = 0.3; = 0.7 model, and the Q = 0.3; = 0.7 model more size evolution than 
the Q = 0.1 model because of the steeper dependence of 1/H[z) at the typical formation 
redshifts of disks to 2; ~ 1. The predictions of both the Prantzos & Aubert (1995) and 
Chiappini et al. (1997) models are quite similar, encouragingly enough given that Prantzos 
& Aubert (1995) and Chiappini et al. (1997) models represent completely independent 
efforts to model the evolution of the Milky Way galaxy. 

Taking R/Vc as the measure of size for a given mass halo, we plot both the Courteau 



- 16 - 



sample ai z = and the higher-redshift data of Vogt et al. (1996,1997) scaled appropriately 
so that the Courteau data is centered on unevolved scale size of a galaxy at 2; = 0. At face 
value, a comparison of the Vogt data with the Courteau data indicates that there has been 
size evolution from z = to z = 1 as any of the models here would predict (Mao et al. 
1998). Nevertheless, the lack of a strong trend in redshift across the Vogt data-set makes 
one suspicious that there may be selection effects at work or even systematic errors in the 
measurements of parameters which may produce the observed differences. In any case, 
strong conclusions must await the compilation of a larger high-redshift dataset, where the 
selection effects have been more carefully quantified. 

We also present model scalings of the number, luminosity, and surface brightness 
expected for L^, galaxies in Figure 1. By construction, the hierarchical models produce more 
number evolution than the infall models, which involve no evolution in number. Except for 
the infall model based upon the Chiappini et al. (1997) prescription, our "Infall" model 
produces similar evolution in luminosity as the hierarchical models to 2; ~ 1. The infall 
models also produce less surface brightness evolution than the hierarchical models. 

4.2. Tully-Fisher Relationship 

To assess hierarchical and infall models, we compare their predicted Tully-Fisher 
relationships with both low and high-redshift observations in the left panels of Figures 2-3. 
For the Monte-Carlo simulations, we use the same selection effects as already specified for 
the low (Pierce & Tully 1988) and high-redshift (Vogt et al. 1996, 1997) samples. We have 
added the Pierce & Tully (1992) fit to these plots for comparison. Naturally, for both our 
hierarchical and infall models, we obtained good agreement with the Pierce & Tully (1988) 
sample since we used that sample to adjust the mass-to-light ratio (we assumed that a 
Mfjj = —21 galaxy had a mass 1.2 ■ IO^^Mq) and its assumed log-normal scatter (one sigma 
scatter of 0.3 dex). At higher redshift, we again obtain basic agreement with the Vogt et al. 
(1996,1997) sample for both our infall and hierarchical models. 

4.3. Size-Vc Relationship 

We also compare our model predictions with the low and high-redshift observations for 
the size-Vc relationship in the right panels of Figures 2-3. Again, we apply the selection 
criteria given in §3 to the low (Courteau 1997) and high-redshift (Vogt et al. 1996,1997) 
model results. At low redshift, we obtain basic agreement with the observations of Courteau 
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(1997) though there seems to be a slight shift in our models toward larger sizes. No 
adjustment of our models was made to obtain agreement with the Courteau (1997) sample, 
and so this comparison can be considered a self-consistency check on our z = models. 

At high redshift, sizes for both models are consistent, if not a little larger than the 
sizes in the Vogt et al. (1996,1997) sample. One of the most surprising thing about a 
comparison of the models with the observations is the significant size evolution observed in 
the lowest redshift bin [z < 0.6) relative to the models. In fact, as discussed in relation 
to Figure 1, the low redshift {z < 0.6) points seem to have undergone more size evolution 
than the high redshift {z > 0.6) points. As the low redshift points are primarily low 
luminosity galaxies and the high redshift high luminosity galaxies, this could point to some 
luminosity-dependent evolutionary trend though the numbers are still too small to make 
any claims toward this end. 



4.4. Size-Magnitude Relationship 

As so often in making low-to-high redshift comparisons, freedom in the choice of 
the z = no-evolution model can be very important in interpreting the high-redshift 
results. In particular, due to the fact that each redshift bin in the Simard et al. (1999) 
sample contains galaxies of a particular luminosity, significant evolution in disk surface 
brightness would appear to be present, simply as a result of correlations between luminosity 
and surface brightness at 2 = 0. There are also important surface brightness selection 
effects in constructing the Simard et al. (1999) sample. We illustrate the importance of 
these considerations in Figure 4 by presenting a no-evolution model, hereafter referred to 
as our fiducial no-evolution model, identical to our z = hierarchical and infall models 
except there is no-evolution in the disk size, number, surface brightness, or luminosity, 
i.e., Ei{z) = Ej^{z) = Eji{z) = 1. We also present the surface brightness distributions 
recovered from a similar no-evolution model differing only in its use of the M^j = —21 
surface brightness distribution for all luminosities. We also present the surface brightness 
distributions recovered both by including a less conservative surface brightness selection 

[ 1, /8i4 + 51ogrrf(//) < 21, 

Supihu,rdi")) = I l-|(/8i4 + 51ogrrf(//)), 21< /814 + 5 logr,(//) < 24, (39) 
i 0, 24< /8i4 + 51ogr,(//) 

(more resembling the one used by Simard et al. (1999)) and without including surface 
brightness (Sup) selection at all. Notice the apparent increase in surface brightness for 
our fiducial model, where SB oc L^/^, as a function of z relative to our constant surface 
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brightness model. It is interesting to note that there is an absence of galaxies at surface 
brightnesses close to the threshold for detection in the Simard et al. (1999) sample. 

With these caveats in mind, we compare our model i?-band surface brightness 
distributions with both the low and high redshift observations in Figure 5. The models 
seem to be in rough agreement with the surface brightness distribution of the observations. 
Given that model populations increase in 5-band surface brightness by ~ 1.5™ to z ~ 1 
(see Figure 1), this suggests a similar increase in the surface brightness of disks to z ~ 1. 
The models themselves show no large differences. As so often, the uncertainties in the z = 
modeling are large enough to preclude detailed discrimination among models, especially 
given the limited high redshift data sets. 

In Figures 6-7, we plot the observed luminosity-size distributions at low (de Jong & 
van der Kruit 1994) and high (Simard et al. 1999) redshift and compare them with those 
obtained for the Q = 0.3/Q\ = 0.7 hierarchical and infall models. We show a similar 
comparison of the hierarchical model with the high redshift Lilly et al. (1998) large disk 
data-set in Figure 8. We also plot the cumulative size and luminosity distributions along 
the vertical and horizontal axes, respectively, for both the observations (histogram) and the 
models (lines). 

As in Figures 4-5, the models predict too many low surface brightness galaxies relative 
to the observations. This also results in too many large model galaxies and too many low 
luminosity model galaxies relative to the observed size and luminosity distributions at 
intermediate to low luminosities, especially for the Lilly et al. (1998) large disk sample even 
with our more conservative surface brightness limit. For large, luminous galaxies, there is 
no obvious change in numbers to high redshift. 

In Figures 6-7, there are also relatively large differences in normalization between the 
observations and models. This is apparently the result of large-scale structure (there are 
small groups/clusters at 2 ~ 0.8 and z ~ 1.0 in the Groth Strip). It is therefore difficult 
to make extremely quantitative statements about the evolution in the number density 
of galaxies with specific surface brightnesses, sizes, and luminosities across the redshift 
intervals surveyed. 

Nevertheless, the abundance of high surface brightness disk galaxies at high redshifts 
relative to the model predictions is surely conspicuous and suggests that there has been a 
significant increase in the number of high surface brightness disk galaxies to z ~ 1 as we 
have already argued in comparing the model and observed surface brightness distributions. 
Again, shifting the surface brightness distribution toward these high surface brightnesses 
(here by the ~ 1.5™ predicted by the models) is an obvious way of accommodating this 
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increase. Of course, any argument based on the normalization of specific galaxy populations 
is subject to considerable uncertainties important when such small contiguous areas are 
being probed. 

To provide a visual comparison between our no-evolution and evolutionary models, we 
include a simulation of a patch of the HDF {Ipsuw: -8^4501^; and Vpeoew) in Figure 9 using 
our fiducial no-evolution model (panel a), the f2 = 1 hierarchical model (panel b), and the 
Q = 0.3,Qa = 0.7 infall model (panel c) for comparison with the actual HDF North and 
South (panel d). Clearly, the lack of high surface brightness galaxies is apparent in the 
no-evolution model relative to the HDF and even somewhat in the apparent in the f2 = 1 
hierarchical model relative to the HDF. Of course, our simulations do not include ellipticals 
or peculiars, so the actual HDF will include more bright objects than the simulations. 

5. Discussion 

There is a real question about a lack of low surface brightness galaxies relative to our 
predictions, especially as compared to the no-evolution model predictions. This conclusion is 
somewhat dependent on the assumed correlation between surface brightness and luminosity 
as is evident in Figure 4. This conclusion is also dependent on the selection biases against 
low surface brightness galaxies not being stronger than those considered here. 

There is an extensive literature discussing surface brightness selection biases (Disney 
1976; Allen & Shu 1979) and various attempts to derive the bivariate luminosity-surface 
brightness distribution of galaxies (McGaugh 1996; Dalcanton et al. 1997b; Sprayberry 
et al. 1997). Surface brightness has a particularly strong effect on isophotal magnitude 
determinations, especially for low surface brightness galaxies; and this can introduce 
significant errors in the magnitude determinations, so the effective volume probed for these 
galaxies is significantly smaller than it is for equivalent luminosity high surface brightness 
galaxies (McGaugh 1996). 

Simard et al. (1999) in a detailed quantification of the selection effects of the DEEP 
sample do not consider the effect of surface brightness on the magnitudes and sizes recovered 
since typical errors were found to be 0.2™ (Simard 1999, private communication). Despite 
the relatively small size of this error, it is not entirely clear to the present authors that the 
errors would not become quite significant for the lowest surface brightness galaxies in the 
sample, particularly those just marginally detectable given the chosen object identification 
and photometric parameters. Secondly, Simard et al. (1999) considers disk galaxies to be 
optically thin whereas the observations of Lilly et al. (1998) are more consistent with disks 
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being optically thick. Highly-inclined optically thin disks would be much more detectable 
than face-on or optically-thick disks. The upshot is that at many apparent magnitudes and 
radii, Simard et al. (1998) would suppose that at least some highly inclined galaxies would 
be detectable and therefore the selection function Sup there would be non-zero when in 
reality if disks were optically thick it would be zero. For these reasons, we used a slightly 
more conservative selection function in surface brightness than that given in Figure 4 of 
Simard et al. (1999) (see §3.2). 

Another possibility, not considered here, is that low surface brightness galaxies might 
form relatively late, meaning that their mass-to-light ratios remain relatively large until 
relatively recent epochs. Of course, prima facie, this would seem unlikely given the 
apparently constant slope in the Tully-Fisher relationship to faint magnitudes. 

In their own analysis of their sample of ~ 200 galaxies, Simard et al. (1999) concluded 
that there had been little evolution in the surface brightness distribution of disk galaxies 
when all selection effects had been carefully considered, quite in contrast to our estimated 
~ 1.5™ of i?-band surface brightness evolution. Little consideration, however, was paid to 
the evolution in the total numbers of high surface brightness galaxies. Here, we find that 
the number of high surface brightness galaxies dramatically exceeds that predicted by the 
evolutionary models considered here, and we have argued that this provides evidence for an 
evolution in the surface brightness distribution of disk galaxies. 

Our interpretation seems to be furthermore supported by the lack of low surface 
brightness galaxies relative to our models. For no-evolution in the disk surface brightness 
distribution really to be present as Simard et al. (1999) claims, high-redshift intervals should 
have similar numbers of low surface brightness galaxies to those found in local samples, and 
these galaxies seem to be deficient, even with respect to our models which show significant 
evolution in surface brightness. 

While the conclusions of Simard et al. (1999) appear to have been carefully drawn, 
we would like to suggest that there are significant uncertainties in their determination of 
the mean surface brightnesses in the lowest redshift intervals and therefore the inferred 
evolution in surface brightness due to the small size of the low redshift samples considered. 
By applying the selection effects from the high-redshift bin identically to all redshift 
intervals, Simard et al. (1999) restricted their analysis to that fraction of disk galaxies 
exceeding the high-redshift surface brightness detection limit. Applying these selection 
criteria uniformly to all low redshift intervals severely pares down the low-redshift samples 
and significantly increases the uncertainty of their average surface brightness measure. 
Given the observed range in observed surface brightness (~ 2mag/arcsec^) and typical 
numbers (~ 5 — 6) for the lowest redshift bins, there is a non-negligible uncertainty in the 
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average surface brightness at low redshift, ~ 0.6™. 

Our estimates of ~ 1.5*" of 5-band surface brightness evolution are somewhat larger 
than that inferred by most authors. Roche et al. (1998) found 0.9"^ of surface brightness 
evolution from z ~ 0.2 to 2; ~ 0.9, Lilly et al. (1998) found 0.8"^ of surface brightness 
evolution in their large disk sample, and Schade et al. (1995,1996a) inferred 1.2™ and 
1.5™ respectively to 2; ~ 0.8. Despite different differential measures of surface brightness 
evolution, most of these samples give similar values for the mean disk surface brightness near 
2 ~ 1: 20.79 ± 0.17 for the Roche et al. (1998) sample (0.65 < z), 19.9 ± 0.2 for the Simard 
et al. (1999) sample (0.9 < 2 < 1.1), 20.7 ± 0.25 for the Lilly et al. (1998) large disk sample 
(0.5 <z < 0.75), 20.2 ± 0.25 for the Schade et al. (1995) sample (0.5 < 2), and 19.8 ± 0.1 
for the Schade et al. (1996) sample (0.5 < z < 1.1). Consequently, differences in the 
surface brightness evolution inferred derive from differences in the z = surface brightness 
distributions assumed. We assume a distribution consistent with the local data of de Jong 
& van der Kruit (1994) and Mathewson et al. (1992) as a baseline for measuring evolution 
with a surface brightness peaking faintward of Freeman's Law (~ 21.75 mag/arcsec^; 
Freeman 1970) while Schade et al. (1995,1996) simply makes reference to Freeman's Law 
(~ 21.65 6jmag/arcsec^). Simard et al. (1999), Lilly et al. (1998), and Roche et al. (1998) 
measure surface brightness evolution differentially across their samples. Possible problems 
here are surface brightness selection effects and limited low-redshift samples. 

6. Summary 

In the present paper, we presented models based on two different approaches for 
predicting the evolution in disk properties: a hierarchical forwards approach, where the 
evolution in disk properties follows from corresponding changes in halo properties, and a 
backwards approach, where the evolution in disk properties follows from an infall model 
providing a close fit to numerous observables for the Milky Way. We normalized the models 
to the local z = observations, we made high-redshift predictions for the models, and we 
compared these predictions with high-redshift observations. 

Our findings are as follows: 

• The hierarchical and infall models predict relatively similar amounts of evolution in 
global properties (size, surface brightness, and luminosity) for disk galaxies to z ~ 1. 
Clearly, discriminating between the models will require a careful look at evolution 
in number (and therefore surveys over a much larger area) and/or measurements of 
certain internal properties, like color, star formation, or metallicity gradients of high 
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redshift disks. 

• There is an apparent lack of low surface brightness galaxies in the high-redshift 
observations of Simard et al. (1999) and Lilly et al. (1998) as compared to model 
predictions based on local observations (Mathewson et al. 1992; de Jong & van der 
Kruit 1994). 

• Our model surface brightness distributions produce relatively good agreement with 
the observations, suggesting that the 5-band surface brightness has evolved by ~ 1.5™" 
from z = to 2; ~ 1 similar to that found in the models. This finding is supported by 
the fact that there is a significantly larger number of high surface brightness galaxies 
than in our model predictions, suggesting that there has been a significant evolution 
in number, most easily accommodated by shifting the mean surface brightness of the 
disk population to higher surface brightnesses. This is contrary to the conclusion 
reached by Simard et al. (1999) based on the same data. 

Here the hierarchical and infall models were presented as competing models to describe 
the evolution in the properties of disk galaxies. If the hierarchical structural paradigm is 
roughly correct as is generally supposed, the infall model simply provides a modification of 
the basic hierarchical scalings to account for the fact that the gas infall rate or star formation 
efficiency is not the same at all radii. In this sort of scenario, if there is an appreciable 
formation of structure at low redshift, a consideration of hierarchical scalings is probably 
more appropriate and if there is not, a consideration of scalings following from infall models 
is probably more appropriate. Obviously, at the present time, a complete incorporation of 
radially dependent gas infall and star formation scenarios into a hierarchical paradigm is 
not merited given the lack of high-redshift data needed to constrain such hybrid models. 
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Fig. 1. — Evolution of size, rest-frame B luminosity, number, and surface brightness for a 
fiducial (M ~ 1.2 ■ lO^^M©) galaxy in our Q = 0.3/Q\ = 0.7 hierarchical model, our 
f2 = 0.1 hierarchical model (long dashed line), our f2 = 1 hierarchical model (dotted line), 
our infall model (dashed line), and the Chiappini et al. (1997) infall model (dotted-dashed 
line). Using the Rd/Vc as a measure of size for a given mass halo, we have added the z = 
Courteau (1997) sample and higher- redshift Vogt et al. (1996, 1997) sample to this plot, our 
scaling the Rd/Vc values so that the Courteau sample had a fiducial scale length of unity at 
z = 0. Both the data and models are presented using Q = 0.3, = 0.7, Hq = 70km/s/Mpc, 
and corrected to unattenuated magnitudes (Tully & Fouque 1985). 
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Fig. 2. — Comparison of the observed Tully-Fisher and size-circular velocity relationships 
as a function of redshifts against the presented Q = 0.3, = 0.7 hierarchical model at low 
and high redshift. The low-redshift Tully-Fisher data is from Pierce & TuUy (1988), the 
high-redshift data is from Vogt et al. (1996, 1997), the superimposed line is the Pierce & 
Tully (1992) Tully-Fisher relationship, and the low-redshift Vc-size data is from Courteau 
(1997). 
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Fig. 3. — Comparison of the observed TuUy-Fisher and radius-circular velocity relationships 
against the presented Q = 0.3, I^a = 0.7 infall model at low and high redshift. The data is 
as in Figure 2. 
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Fig. 4. — Comparison of the observed rest-frame B surface brightness distributions 
(histogram) with those from our fiducial no-evolution model (solid line), our fiducial no- 
evolution model with constant Mbj = —21 surface brightness distribution (dotted line), 
our fiducial no-evolution model without surface brightness selection (dashed line), and our 
fiducial no-evolution model with the less conservative selection function of Simard et al. 
(1999) (long dashed line). This illustrates the possible importance of surface brightness 
selection and an assumed luminosity-surface brightness correlation on the conclusions 
derived. The local [z < 0.1) data is from de Jong & van der Kruit (1994) and the high- 
redshift {z > 0.1) data is from Simard et al. (1999). 
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Fig. 5. — Comparison of the observed rest-frame S-band surface brightness distributions 
(histogram) with those from our hierarchical models {Q = 0.3, Qa = 0.7/solid line; 
Q = 0.1/dashed line; Q = 1/dotted line), our infall models (thick dotted line), and our 
fiducial no-evolution model (thick solid line) presented here. The local {z < 0.1) data is 
from de Jong & van der Kruit (1994) and the high-redshift {z > 0.1) data is from Simard et 
al. (1999). 
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Fig. 6. — Comparison of the observed magnitude-radius (rest-frame B) distributions with 
the hierarchical models presented here. The low-redshift data (filled circles) is from de 
Jong & van der Kruit (1994) and the high-redshift data (filled circles) is from Simard et 
al. (1999), and the small dots are the results for the Q = 0.3, Q\ = 0.7 hierarchical model. 
Cumulative size and luminosity distributions are presented on the vertical and horizontal 
axes, respectively, for our = 0.1 hierarchical model (dashed line), our Q = 1 hierarchical 
model (dotted line), our Q = 0.3,Q\ = 0.7 hierarchical model (solid line), and our fiducial no- 
evolution model (thick solid line) for comparison with the observations (histogram). Both 
the data and models are presented using Q = 0.3, fl^ = 0.7, Hq = 70km/s/Mpc, and 
assuming an inclination of 70 deg (TuUy & Fouque 1985). 
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Fig. 7. — Comparison of the observed magnitude-radius (rest-frame B) distributions with 
the Q = 0.3, Qa = 0.7 infall model presented here (§2.2). The data is as on Fig 5, and the 
small dots trace out the model distribution. Cumulative size and luminosity distributions 
are presented for the infall model (solid lines), the no-evolution models (thick line), and the 
observations (histogram) . 
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Fig. 8. — Comparison of the observed magnitude-radius (rest-frame B) distributions from 
the Lilly et al. (1998) large disk sample with the hierarchical models presented here. The 
models and data are as in Figure 6. Both the data and models are presented using Q = 0.3, 
= 0.7, Hq = 70km/s/Mpc, and assuming an inclination of 70 deg (Tully & Fouque 
1985). There are an excess of model galaxies at low magnitudes and large sizes relative to 
the observations. 



Fig. 9. — Panels (a), (b), (c), and (d) show 60" x 80" color images {Ipsuw, Bp^^ow, 
and VpQOQw) for an HDF-depth simulation using our fiducial no-evolution model, an HDF- 
depth simulation using our f2 = 1 hierarchical model, an HDF-depth simulation using our 
Q = 0.3,flx = 0.7 infall model, and a portion of the actual HDF North and South. Clearly, 
our no-evolution model has fewer high surface brightness galaxies than the HDF. Our f2 = 1 
hierarchical model also appears to lack high surface brightness galaxies though the fact that 
we did not include peculiars and bright ellipticals in our simulations would tend to bias the 
eye. 



